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Abstract 



We discuss the effect of electromagnetic environment on laboratory mea- 
surements of the nuclear fusion reactions of astrophysical interest. The radi- 
ation field is eliminated using the path integral formalism in order to obtain 
the influence functional, which we evaluate in the semi-classical approxima- 
tion. We show that enhancement of the tunneling probability due to the 
radiation correction is extremely small and does not resolve the longstanding 
problem that the observed electron screening effect is significantly larger than 
theoretical predictions. 
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Nuclear fusion reactions measured at a laboratory at very low incident energies are 
subjected to the electron screening effect, which originates from bound electrons in the 
target atom (or molecule). The electrons shield the Coulomb potential between the projectile 
and the target nuclei, and it is expected that fusion cross sections are enhanced over the 
prediction which does not take into account this correction [[!]]. Indeed, since the seminal 
work for the 3 He(d,p) 4 He reaction ||, the experimental fusion cross sections (or equivalently 
the astrophysical S factors) are shown to be enhanced as compared to those extrapolated 
from the high energy region where the electron screening effect is negligible down to the low 
energy regime (see Ref. |3[] for a recent review). 

In the low energy limit, the adiabatic approximation is validated, and the electron screen- 
ing effect can be expressed as a constant energy shift of the Coulomb potential, where the 
amount of the shift is given by a difference of electron binding energies between the uni- 
fied and the isolated systems [|lj. The experimental data, however, systematically indicate 
that a significantly larger energy shift is required in order to account for the low energy 
enhancement of the fusion cross sections [|J. This has been a big surprise and also a rather 
puzzling situation, since the adiabatic approximation should provide the upper limit of the 
tunneling probability and thus the upper limit of the energy shift 0-0. No satisfactory 
explanation has been proposed to reconcile this problem so far. Furthermore, the recent 
careful measurement for the stopping power P|| as well as the recent attempts which used 
the Trojan- horse method flnj to measure the bare cross sections at low energies [11, 12 1 have 
re-confirmed that the enhancement of fusion cross section is much larger than the model 
calculation. 

Balantekin et al. studied several effects beyond the electron screening on astrophysical 
reactions, which include vacuum polarization, relativity, bremsstrahlung, and atomic po- 
larization [113. They found that all of these effects are too small to explain the difference 



between the measured and the predicted electron screening energies. Obviously the low en- 
ergy enhancement of the astrophysical reactions has a different origin. Recently, Flambaum 
and Zelevinsky suggested that the virtual photon emission during tunneling may increase the 



penetrability |]14 |. Within a few approximations, which include the closure approximation 



for the relative motion between the projectile and the target nuclei, they derived a static po- 
tential shift due to the radiation correction which is proportional to the second derivative of 
the bare internucleus potential. The second derivative is exactly zero for the point Coulomb 
potential and the potential renormalization which Flambaum and Zelevinsky derived has a 
finite value only well inside the Coulomb potential. Note that such potential renormalization 
can be well compensated with a choice of a nuclear potential between the projectile and the 
target unless the energy dependence is very strong, as in the case of coupling to high-lying 
states in heavy-ion subbarrier fusion reactions fL5l , and it is thus not easy to separate out 
the radiation effect in a clear way. Since the closure approximation may be too crude for 
astrophysical reactions, it is necessary to re-examine the effect of virtual photon emission 
without resorting to the approximations which Flambaum and Zelevinsky used in order to 
assess the role of radiation field in realistic systems. 

The purpose of this paper is to study the effect of coupling between the tunneling motion 
and the electromagnetic field taking a different approach from Flambaum and Zelevinsky. 
To this end, we employ the path-integral formalism for multi-dimensional tunneling [16| 



and evaluate the influence functional in the semi-classical approximation. Besides the semi- 



2 



classical and the dipole approximations, our formalism is exact, and the effect is finite even 
with the pure Coulomb potential. The path integral approach allows us to discuss both 
virtual and real photon emissions during tunneling on the equal footing. In this respect, our 
study is relevant also to the bremsstrahlung in a decay, which has attracted much interest 



for the past few years [[17] [H 

Consider the tunneling motion in the presence of the radiation field. The classical La- 
grangian for the system reads pp|pT | 
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where a represents the polarization index and ej^ is the corresponding polarization vector. 
We employ the Coulomb gauge, and thus the polarization vectors are orthogonal to the 
photon momentum k. We have used the dipole approximation, and introduced the photon 
coordinate Q ka {t), which is independent of the relative motion r between the projectile and 
the target nuclei, egi is the El effective charge given by (ApZx — AtZp) /{Ap + At), At 
and Zt being the mass and the atomic numbers for the target nucleus, respectively, and 
similar for the projectile also. \i is the reduced mass for the relative motion r, and V(r) is 
the potential between the projectile and the target in the absence of the radiation field. As 
usual, we have neglected a term which is proportional to e| x in the classical Lagrangian. 

Our interest is to compute the transition amplitude for the relative motion from r« to 
rf while for the radiation field from the vacuum state |0) to the n-photon state |n) for a 
given incident energy E. This is expressed in terms of the energy representation of the path 
integral amplitude as 



g J™ dT e lET/n J V[r{t)]e (l/h) So dtL °(n J V[Q ka {t)]e (t/h) So dt ^ 
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where the kinematical factor g is proportional to yE/fi fll6| . Here the integral for r is 
carried out for all paths which connect r(0) = r» and r(T) = rf. We square this amplitude 
and sum over n in order to obtain the total transition probability for the relative motion: 

P = Y.\Tn\ 2 (6) 
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where p is the two-time influence functional given by |TJ 



p(r(t),T;r,T) = (0 
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For the Lagrangian given in Eqs. (HD - (@), the two-time influence functional can be expressed 
analytically and is given by 



p(r(t),T;r,T) = exp 
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where f ka and f ka is defined by 
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In the exponent in the influence functional (|Sp, the last term is the real photon contribution, 
while the second and the third terms represent the virtual photon emission p0 |. 

In the radiation problems such as the Lamb shift calculation, it is essential to separate 
out a divergent contribution due to the mass renormalization in order to obtain a physical 
result. For our problem, this can be done by performing a partial integration for the second 
and the third terms in Eq.(P): 
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The first term is proportional to r 2 and is nothing but the mass renormalization, which has 
already included in the kinetic term in Eq. (|2]). The second term vanishes if we choose 
t — at the outer classical turning point (note that we are going to deal with the tunneling 
problem in the following). Retaining only the third term, the regularized influence functional 
p rcg then reads 



p reg (r(t),T;r,T) = exp 
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Here, we have dropped also the — |a;(T — T) term, since it merely changes the definition of 
the energy of the vacuum state. 

We now introduce the semi-classical approximation to the path integrals with respect to 
r as well as to the time integrals in Eq. (|7|) [|,^,^3|]. This leads to a classical tunneling path 
from the outer to the inner turning points along the imaginary time axis. For simplicity, 
in the following, we consider only a head-on collision. Since the coupling strength e^/hc 
is small, we determine the tunneling path by disregarding the radiation field. The classical 
path thus obeys the Newtonian equation with the inverted potential, 
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The influence functional is then evaluated along the classical path. After carrying out 
the angle integral for the photon momentum k, the penetrability (0) in the semi-classical 
approximation thus reads 



where 
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is the regularized two-time influence functional along the classical trajectory z c \. In obtaining 
/, we have introduced the imaginary time evolution, t — > —ir,s — > —it', and T — > —iT , 
and neglected the excitation of the radiation field prior to the tunneling. T is the Euclidean 
tunneling time, where z c i(T ) corresponds to the inner turning point R while z c i(0) is the 
outer turning point Z P Z T e 2 /2E. 

Let us now numerically estimate the enhancement factor f = P/P for the <i+ 3 He reac- 
tion. To this end, we consider the pure Coulomb potential V(r) = ZpZj-e 2 jr from the outer 
turning point ZpZj-e 2 / E to the inner turning point i?o=4.3 fm. The Newtonian equation 
(|j~3l) can then be solved analytically as, 
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Because of the exponential factor, the k integral in Eq. (IT7]) is quickly damped as a function 
of k, and it is not necessary to introduce the momentum cut-off factor for our tunneling 
problem. Figure 1 shows the deviation of the enhancement factor / from unity as a func- 
tion of the incident energy. We see that the enhancement due to the radiation coupling in 
the tunneling region is extremely small and does not play any important role in the astro- 
physical fusion reaction. The dashed line shows the virtual photon contribution separately, 
which is obtained by neglecting the second term in Eq. flT^) . One finds that the virtual 
photon emission provides the dominant contribution, although the real photon emission is 
not negligible. Table 1 shows the enhancement factor for several reactions. We take the 
same value for the inner turning point i?o as in Ref . |13] . In all the reactions considered, the 
effects of the radiation coupling during tunneling is almost negligible and the enhancement 
factor never exceeds 10~ 3 %. This effect is even much smaller than the other small effects 
considered in Ref. [0 . We thus conclude that the radiation correction is not helpful for the 



low energy enhancement of nuclear astrophysical reactions. 
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In summary, we have studied the effect of the radiation coupling on the tunneling motion 
of charged particles. To this end, we have employed the semi-classical approximation to the 
transition amplitude in the path integral representation. We have shown that the effect 
is almost negligible for nuclear astrophysical reactions, and the large electron screening 
problem remains unsolved. The radiation correction which we discussed in this paper is 
even smaller than the vacuum polarization effect, and would be negligible as compared to 
the latter. The origin of the discrepancy between the measured and the calculated electron 
screening energies still seems to be an open problem. 

We thank the European Centre for Theoretical Studies in Nuclear Physics and Related 
Areas (ECT*), where a part of this work was done, for its hospitality and for partial support 
for this project. This work was supported in part by the U.S. National Science Foundation 
Grant No. PHY-0070161 at the University of Wisconsin, and in part by the University 
of Wisconsin Research Committee with funds granted by the Wisconsin Alumni Research 
Foundation. 
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E (keV) 

FIG. 1. The deviation of the enhancement factor / due to the radiation coupling from unity 
as a function of the incident energy for the <i+ 3 He fusion reaction. The solid line includes both the 
real and the virtual photon emissions during tunneling, whereas the dashed line represents only 
the virtual photon contribution. 
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TABLES 

TABLE I. The enhancement factor / for several reactions obtained at the lowest experimental 
energy E m \ n and with the inner turning point Rq indicated. 



Reaction 


E min (keV) 


Ro (fm) 


/-I 


3 He(d,p) 4 He 


5.88 


4.3 


2.744 x 10" 7 


D( 3 He,p) 4 He 


5.38 


4.3 


2.731 x 10~ 7 


6 Li(p,a) 3 He 


10.74 


3.0 


3.758 x 10~ 6 


7 Li(p,a) 4 He 


12.70 


4.3 


3.598 x 10~ 6 


H( 6 Li,a) 3 He 


10.94 


3.0 


3.761 x 10~ 6 


H( 7 Li,a) 4 He 


12.97 


4.3 


3.602 x 10~ 6 


10 B(p,a) 7 Be 


18.70 


3.3 


6.095 x 10~ 6 


n B(p,a) 8 Be 


16.70 


2.0 


1.169 x 10~ 5 
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